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Abstract: 

Coupled discrete models abound in several areas of physics. Here we provide an extensive set of exact 
quasiperiodic solutions of a number of coupled discrete models in terms of Lame polynomials of order one 
and two. Some of the models discussed are (i) coupled Salerno model, (ii) coupled Ablowitz-Ladik model, 
(iii) coupled saturated discrete nonlinear Schrodinger equation, (iv) coupled 0^ model, and (v) coupled 4)^ 
model. Furthermore, we show that most of these coupled models in fact also possess an even broader class 
of exact solutions. 
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1 Introduction 

In a recent paper fT] we have obtained solutions of a coupled (p^ and and a coupled 0^ model in terms 
of Lame polynomials of order one and two even though most of these solutions are not the solutions of 
the corresponding uncoupled problem. The purpose of the present paper is to carry out a similar study 
for a number of coupled discrete field theory models. In particular, we obtain exact solutions of a (i) 
coupled Ablowitz-Ladik (AL) model, (ii) coupled Salerno model, (iii) coupled saturated discrete nonlinear 
Schrodinger equation (DNLSE), (iv) coupled (j)^ model, and (v) coupled 0^ model. We also show that 
unlike the continuum field theory models, many of the discrete coupled field theory models possess an even 
broader class of exact solutions. Moreover, we show that as in the uncoupled case [2], even the coupled 
AL, coupled Salerno and coupled DNLSE models follow from the same Hamiltonian but with a different 
Poisson bracket (PB) structure. 

The motivation for this work comes from the fact that there are many physical situations where a dis- 
crete field theory is appropriate to model the phenomena of interest with a specific coupling between the 
two fields. One such phenomenon of current intense interest is the coexistence of magnetism and ferroelec- 
tricity (i.e. magnetoelectricity) in a given material. This is a highly desired functionality in technological 
applications involving cross-field response, switching and actuation. In general, this phenomenon is re- 
ferred to as multiferroic behavior 0. Recently, two different classes of (single phase) multiferroics, namely 
the orthorhombically distorted perovskites [1] and rare earth hexagonal structures [5], have emerged. The 
latter show magnetic domain walls in the basal planes which can be modeled by a coupled 0^ model [6] in 
the presence of a magnetic field. Coupled 0^ models 13 [HI [9] also arise in the context of many ferroelectric 
and other second order phase transitions. The coupled (/)^ model for multiferroics [6] has a biquadratic 
coupling whereas the coupled 0^ model for a surface phase transition with hydration forces [9] , relevant in 
biophysics context, has a bilinear coupling. Other types of couplings are also known for structural phase 
transitions with strain |10j . 

Examples of coupled discrete AL, coupled discrete Salerno and coupled saturated DNLS models are 
also known [Til HI US]- Similarly, there are analogous coupled models in field theory [HI [15]. Several 
related models have been discussed in the literature and their soliton solutions have been found \16\ [T7\ 
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[H dg [201 [211 [221 [23] including periodic ones [21 [23 [26] . 

The paper is organized as follows. In Sec. II we first show that the coupled AL, coupled Salerno and 
coupled DNLSE models can all be obtained from the same Hamiltonian but with different PB structure. 
We also obtain additional conserved quantities in these models. In Sec. Ill we provide the solutions for 
the coupled Salerno model in terms of Lame polynomials of order one and two as well as a broader class 
of solutions. In Sec. IV we provide similar solutions for the coupled AL model. In Sec. V we show that 
unlike the coupled Salerno and coupled AL cases, the coupled saturated DNLSE model only admits Lame 
polynomial [27] solutions of order one but not of order two. Besides, we have not been able to obtain a 
broader class of solutions in this case. Section VI is devoted to solutions of a coupled discrete <f)^ model 
[Tll28j in terms of Lame polynomials of order one and two and also a broader class of solutions. In Sec. VII 
we discuss solutions of a coupled (j)'^ model introduced by us recently [22] in terms of Lame polynomials 
of order two as well as a broader class of solutions. Note that the solutions in terms of Lame polynomials 
of order one have already been obtained by us in |29j . Section VIII contains the summary of main results 
and possible future directions. 

2 The Model for Coupled saturated DNLSE, coupled AL and coupled 
Salerno Equations 

We have previously shown [2] that the uncoupled Salerno model [30], the uncoupled AL model [31j and 
the uncoupled saturated DNLSE model can all be deduced from the same Hamiltonian 

N 



\Un-Un+l\ \Un\ + ^ ln(l + ^ n„ 



(1) 



n=l 

but with different PB structure. We now show that even the coupled Salerno, coupled AL and the coupled 
saturated DNLSE models can all be derived from the same Hamiltonian given by 

N r 



n=l 



\Un - n„+ip + \Vn - Vn+lP " — l^nP 



Z^2 I 1 2 ^1 T / 1 2 l2^ 

\Vn\ + -n ln(l + )Ul + lJ-2\Vn\ J 

1^2 fJ-i 



(2) 
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with the equations of motion in the two field variables Un and in all three cases being 



iUn = [Un, H] , iVn = [Vn, H] 



(3) 



The difference in the equations of motion comes from a different definition of the PB and consequently a 
different definition of the time derivative. The PB structure in all three cases can be compactly written as 



N r 



[u, v] = y: 



n=l 



dU dV dU dV dU dV dU dV 

+ 



dun 9u* 9m* dun dvn dv^ dv^ dvn 



[I + XlluJ^ + X2\Vnn ■ 



(4) 



Coupled saturated DNLSE 

On using Eqs. ([2]) to (jl]) with Ai = A2 = yields the coupled saturated DNLS equations 



idUn/dt + [Un+l + Un-l - 2li„] + 



n I jUn 
+ /il|UnP + H2\Vn\'^) 



0. 



idVn/dt + [Vn+l + Vn-1 - 2Vn] + 



[^2 - ^) Vn + J^2(A*1 I Un ^ +M2 | Vn 



0. 



(5) 



(6) 



/^2(1 + Wl^nP + f^2\Vn\'^) 

It is easily checked that in this case, apart from the Hamiltonian ([1]), two other conserved quantities are 
power Pu and Py defined by 



N 



N 



Pu — ^ ^ \Un\ ) Pv — ^ ^ \Vn\ 



(7) 



n=l 



n=l 



Coupled Salerno Model 

If instead, we use Eqs. ([2]) to (HD with Ai = ni and A2 = /i2 then we obtain the coupled Salerno model 
with field equations 



idUn/dt + [m„+i + Un-l - 2n„] + (/Ui | U„ ^ +^2 I 1^) 



z/1 — 

-Un+l + tt„_l H Ur, 



0, (8) 



idvn/dt + 



t;„+i + Vn-i - 2 + 



+ (pi I Un P +//2 I t-n 1^) 



1^2 — 2/i2 

+ Vn-1 -\ Vr, 



= 0. 

(9) 



It is easily checked that in this case, apart from the Hamiltonian ([T]), the other conserved quantity is power 
P given by 

N 

P=Y,W + f^lK\'^+f^2\Vn\^]. (10) 

Coupled AL Model 



n=l 
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In the special case when vi = 2/ii and f2 = 2/^2, then the coupled Salerno model reduces to the coupled 
AL model with the field equations 

idUn/dt + [Un+l + Un-1 - 2n„] + (^1 I Un ^ +^2 I Vn + U^-l] = , (11) 

+ (/Wl I Un P +M2 I Vn |^)K+1 + Vn-l] = 0. (12) 

It is interesting to note that in this case, apart from the Hamiltonian ([T]) and power P as given by Eq. 
()10p . generalized momentum Pm given by 

N 

Pm=^ ■i[fJ'l{UnUn+l - ■U*n„+i) + fl2{VnVn+l " K^n+l)] , (13) 
n=l 

is also conserved. 

One remark is in order here. Just as the uncoupled Salerno model interpolates between AL and DNLSE, 
it is easy to see that the coupled Salerno model as given by Eqs. ([5]) and Q also interpolates between 
coupled AL model (as given by Eqs. (|11|) and (|12p ) and coupled DNLSE. In particular, in the limit i^i = 2/ii 
and 1^2 = 2/U2, the coupled Salerno model Eqs. ([5]) and Q go over to the coupled AL model Eqs. ([TT]) 
and (|12p . On the other hand, in the limit /^i = /i2 = but with H2/lJ'i = c and z^2/^i = '^i^i coupled 
Salerno model Eqs. ([SD and @ go over to the coupled DNLS equations 

idUn/dt + [Un+l + Un-l - 2n„] + l^id Un ^ +C | Vn P)n„ = , (14) 
idVn/dt + [Vn+1 + Vn-l " 2Un] + Z^ic(| ti„ ^ +C | V„ P)v„ = . (15) 

3 Solutions of the Coupled Salerno Model 

We now show that the coupled Salerno model as given by Eqs. ([HD and ^ has Lame polynomial solutions 
of order one as well as of order two. In fact, it turns out that the coupled model has an even broader 
class of exact solutions of which Lame polynomial solutions of order one and two are the special cases. We 
remind that so far as we are aware of, the uncoupled Salerno model has no known exact solutions. 
We start with the ansatz 

n„ = /ncxp [-«(a;it + (5i)] , Vn = gneyip[-i{uj2t + 62)] , (16) 
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idvn/dt + 



2^2 

Vn+1 + Vn-l Vn 

Ml 



with /„ and gn satisfying 

f^ + agl = b, a,b>0. (17) 

Here a, 6 > are two positive numbers while 5i, 62 are two arbitrary parameters. On substituting this 
ansatz in Eqs. ([8]) and ([9]) we find that this is a consistent ansatz provided 

a= — , b = , 0Ji = — , 0J2 = — n- . (18) 

/"I /"I w ^J'l 

This imphes that such solutions are possible only if /ii,/i2 < 0. Further, since 

i^ = ^^. (19) 

which is a real number, hence the solutions obtained from here are in general only quasiperiodic. Only if 
is a rational number, will the solutions be periodic. 

Clearly this is a very general ansatz which admits a broad class of solutions including Lame polynomials 
of order one and two. 

Lame polynomial solutions of order one 

(i) One solution is 

/n = ^dn[/3(n + C2),m] , Qn = B^/msn[P{n + C2),m] , (20) 

provided Eq. ([TH]) is satisfied and further 

b = A^, fiiA^ = fi2B\ (21) 

Note that /3 is completely arbitrary. Using the fact that dn(x, m) has period 2K{m) while cn(x, m) and 
sn(a;,m) are periodic functions with period 4i^(m), it follows that for the solution (|20p. Un,Vn satisfy the 
boundary condition 

, 2K(m) = Un , V AK(jn) = Vn ■ (22) 

Here K{m) is the complete integral of the first kind. 

(ii) Another solution is 

fn = A^/mcn[l3{n + C2),m], Qn = B^/msn[f3{n + C2),m] , (23) 
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provided Eq. p8]) is satisfied and further 

b = mA^ , fiiA"^ = fi2B^ . 
For tlie solution (|23l) . Un,Vn satisfy the boundary condition 

W , 4K(m) = Un , V 4K(rn) = Vji ■ 

In the limit m = 1, both these solutions go over to the hyperbolic solution 

/„ = Asech[l3{n + C2)] , Qn = Btanh[f3{n + C2)] , 

Lame polynomial solutions of order two 

(iii) One solution is 

fn = Adn^[l3{n + C2), m] + B , gn = F^/msn[^{n + C2), m]dn[/3(n + C2),m 



provided Eq. (jlSp is satisfied and further 



a2 

b=—, fxiA^ = fi2F\ A = -2B. 



For the solution (j27|) . Un,Vn satisfy the boundary condition (j22|) . 
(iv) Another solution is 

/„ = Adn'^[l3{n + C2),m] + B , gn = Fmsn[/3(n + C2), m]cn[/3(n + C2), m] 
provided Eq. (fT8]) is satisfied and further 

^^m^^ ^^^2 ^^^^2^ (2-m)^ = -2B. 

For the solution (j29]) . Un,Vn satisfy the boundary condition 

U 2K(m) = Un , , 2K(m) = fn • 

n-l- ^ n-l- ^ 

In the limit m = 1, both solutions (p7|) and ([29|) go over to the hyperbolic solution 
/„ = Asech'^[l3{n + C2)] + B , 5(n = -Ftanh[/3(n + C2)]sech[/3(n + C2)] . 



(v) Apart from these, several other solutions are possible. For example one can have nonperiodic 
solutions like 

f ^ __Bn 

provided Eq. (fT8|) is satisfied and further 

h = A^, ^ilA'^ = ^I2B\ (34) 
One can obviously write down a wider class of such solutions. For example 

provided Eqs. (fTHI) and (fM|) are satisfied. 

(vi) Yet another possible periodic solution is 

fn = Acos[(3{n + C2)], gn = Bsm[f3{n + C2)], (36) 
provided Eq. (fTSl) is satisfied and further 

6 = yl2, ;uiA2 = /i2S2. (37) 
In this case both Un,Vn satisfy the periodicity condition 

= Un , V^^2. = Vn ■ (38) 

It turns out that apart from the general solution as given by Eqs. (jl6p and ()17p . there is another 
possible general solution given by 

Un = /nexp [-^(wit + (5i)] , u„ = 5t„exp [-i(tJ2i + ^2)] , (39) 

but now /„ and gn satisfy 

fl-cigl = h, 0,5 >0. (40) 
On substituting this ansatz in Eqs. ^ and Q we find that this is a consistent ansatz provided 

a = , b= , uji = — , UJ2 = — 2~ • (41) 

m fii m fif 
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This implies that such solutions are possible only if fii, ^2 have opposite signs. Clearly this is a very general 
ansatz which admits a broad class of solutions. As an illustration we discuss a few such solutions, 
(vii) One solution is 

J ^ ^ ^ ^ B^/msn[P{n + C2),m] ^^^^ 

dn[/3(n + C2), m] ' dn[/3(n + C2), m] ' 

provided Eq. ()4ip is satisfied and further 

b = A'^, /ii<0, H2>0, =/i25^ (43) 



In this case both m„, v„ satisfy the periodicity condition (j22p . Note that if we interchange fn and then 
Hi > 0, fi2 < 0. In the limit m = 1, this solution goes over to the hyperbolic solution 

/„ = Acosh[/3(n + C2)], gn = B smh[f3{n + C2)] , (44) 

(viii) Another solution is 

^ ^ ^ ^ Fy/^sn[/3{n + C2),m] ^^^^ 

dn^[/3(n + C2), m] ' dn^[/3(n + C2), m] ' 

provided Eq. (|1T]) is satisfied and further 

6 = — , /ii<0, fi2>0, \fii\A'^ = fi2F^ , A = -2B. (46) 

In this case both Un,Vn satisfy the periodicity condition (j22p . In the limit m = 1, this solution goes over 
to the hyperbolic solution 

/„ = Acosh2[/3(n + C2)] , 5(„ = Bsinh[^(n + C2)]cosh[/3(n + C2)] , (47) 

(ix) Apart from these, several other solutions are possible. For example one can have the following 
nonperiodic solution 

provided Eq. (|1T]) is satisfied and further 

6 = yl^ /ii<0, /X2>0, ^ //2i^2_ (49) 

One can, easily write down a wider class of such solutions. 
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4 Solutions of the Coupled AL Model 

We show that for the coupled AL model characterized by Eqs. (|lip and (jl2p . one not only has solutions 
similar to those in the previous section (for the coupled Salerno case), but just like the uncoupled AL case 
|32| . even coupled AL equations have moving periodic solutions in terms of Lame polynomials of order one. 

As in the previous section, if we start with the ansatz as given by Eqs. p6|) and pzp or Eqs. (|39p and 
(j40p . then it is easy to show that the entire discussion of the previous section goes through except that 
since in the coupled AL model z^i = 2/ii,z^2 = 2/^2, hence in the coupled AL model with the above two 
ansatze, uji = 2, U2 = But for this minor change, all the nine solutions given in the previous section 
are also solutions of the coupled AL model under the identical conditions (except ui = 2,uj2 = 

In addition to these nine solutions, we now show that as in the uncoupled case |32j, the coupled AL 
model also admits moving periodic solutions in terms of Lame polynomials of order one. 

(i) For example, it admits mixed moving periodic kink-pulse solution 

Un = Aexp[—i{u;it — kin + (5i)]dn[/3(n — vt + 62), m] , 

Vn = B exp[—i{uj2t — k2n + 53)]-v/msn[/3(n — vt + 82), m] , (50) 

provided 



wi = 2 
a;2 = 2 



l-(l + /i25 



2, cos(/ci)dn(/3, m) 



cn2(/3, m) 



M , j,2NC°s(A;2)dn(/3,m) 



/ii cn(/3,m) 
1 = ij.iA'^cs^{l3,m) - fi2B'^ns^{l3,m) , 

/?.. = ^=(Mil±if2^, ,„(ft,„, = 2EM, (51) 

cs(p,mj sm(Ki) 

where cs(/3,m) = cn(/3, m)/sn(/3, m) and ns(/3,m) = l/sn(/3,m). Note that since it is a moving periodic 
kink-pulse solution, it must not only satisfy the periodicity condition (i22]l but it must also satisfy the 
periodicity condition 

U^,2n = Un , V^,2n = Vn ■ (52) 

The periodicity conditions (j22p and (j52p imply that Un,Vn are periodic solutions provided there exist 
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integers ni,n2, 77.3,714 such that 



2K(m) 2tt 4K(m) 27r 
ni — =n2— , ns — = ua— . 

(3 ki p k2 

(ii) Another mixed moving periodic kink-pulse solution that it admits is 



(53) 



Un = Aeyi\>[—i{u)it — kin + 5l)]^/rnc\l[j3{n — vt + 82), 'm] , 
Vn = -Bexp[-z(a;2t - k2n + 5-i)\^/rnsn[l5{n - vt + 82), m] , 



(54) 



provided 



uji = 2 



UJ2 = 2 



1 - (1 + m^i2B 



2^ cos(/ci)cn(/3, m) 
dn^(/3, m) 



dn(/3, m) 

1 = ^llA^ds^{(3, m) - fi2B^ns'^i/3, m) , 
2 s'm{ki) 2i 



■[l + m/i25^], dn(^,m) 



sin(A;2 



(55) 



ds(/3,m)'' ' ■' ' ^f^' ^ sin(A;i) ' 

where ds(/3,m) = dn(/3, m)/sn(/3, m). Note that since it is a moving periodic kink-pulse solution, it must 
not only satisfy the periodicity condition (I25p but it must also satisfy the periodicity condition (I52p . The 
periodicity conditions ()25p and (I52p imply that Un,Vn are periodic solutions provided there exist integers 
ni , n2 , na , 77,4 such that 



4i^(m) 



ni- 



/3 



27r 4E:(m) 27r 

= 724 — 
K2 



In the limit m = 1, both these solutions reduce to the moving pulse- kink solution 

Un = A exp[— i(a;it — kin + 5i)]sech.[f3{n — vt + 82)] , 
Vn = B exp[—i{uj2t — k2n + 5s)] tanh[/3(n — vt + 62)] , 



(56) 



(57) 



provided 



sinh2(/3) = i^iA^ - ^2S^cosh2(/3) , 

uji = 2[l - (1 + fi2B^) cos(A;i) cosh(/3)] , a;2 = 2 



^ - (1 + ^2^2^ COs(fe2) 



vf3 = 2(1 + fi2B^) sin(A;i) sinh(/3) , ^^^M = sech(/3) . 

sm(Kij 
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(58) 



Notice that in case ki = k2 = v = 0, the solutions (|50|) . (jMj) and (|57j) become stationary coupled, 
periodic pulse-kink solutions provided relations (|5ip . (|55|) and (|58|) with fci = A;2 = = are satisfied. 
However, we have already shown (and it can also be verified from relations (|5ip . (|55|) and (j58|) ) that the 
solutions (|50|) . (|54|) and (|57|) with /ci = A;2 = ^ = also hold good under the stronger conditions as given 
by Eqs. (USD, (gl]) and ([Ml) with ui = 2^i, 1^2 = 2/X2. 

(iii) It also admits a coupled moving periodic pulse solution 



Un = ^exp[— z(a;it — kin + (5i)]dn[/3(n — vt + 62), m] 
Vn = B exp[-i{u>2t - k2n + 5^)](ln[j3{n - vt + 62), m] 



(59) 



provided 



ki = k2 , oJi = 2 



cos(/ci)dn(/3, m) 



, ^2 = 2 



/i2 cos(A;i)dn(/3, m) 



cn2(/3, m) 



2 sin(A;i) 



(60) 



cs(/3, m) 

Note that since it is a moving periodic pulse solution, it must not only satisfy the periodicity condition 
(j3ip but it must also satisfy the periodicity condition (j52p . The periodicity conditions (j3ip and (|52p imply 
that Un,Vn are periodic solutions provided there exist integers ^1,^2,7^3, 71.4 such that 



2/C(m) 27r 2K(m) 27r 
ni — = ^2— , ns — = n4— 

p fci /? fc2 



(61) 



(iv) Another coupled periodic moving pulse solution is 



Un = Aexp[—i{ujit — kin + 6i)]\'7ncn[f3{n — vt + 62), m] , 
Vn = B exp[— i(a;2t — ^2"- + (^3)]\/mcn[/3(n — vt -\- 82), ni] , 



(62) 



provided 



ki = k2 , uji = 2 



1 



cos(/ci)cn(/3, m) 



dn {I3,m) 



UJ2 



/i2 cos(A;i)cn(/3, m) 
/ii dn^(/3,m) 



1 = (^i^2^^252)ds2(/3,m), /3t; 



2 sin(A;i 



(63) 



ds(/3, m) 

Note that since it is a moving periodic pulse solution, it must not only satisfy the periodicity condition 
(|25p but it must also satisfy the periodicity condition (j52p . The periodicity conditions ()25p and (j52p imply 
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that Un,Vn are periodic solutions provided there exist integers rai, 712,^3, 714 such that the condition (|56|) is 
satisfied. 

(v) Finally, it also admits a mixed coupled moving periodic pulse solution 



Un = ^exp[— ^(cjit — kin + 5i)]dn[/3(n — vt + 82), m] , 
Vn = B eyi.Y>[—i{uj2t — + 53)]y/rncn[(3{n — vt + 82), m] 



(64) 



provided 



a;i-2 = - 1-(1- m)^2-B^ 



2 cos{ki)dn{/3 , m) 



cn^(/3, m) 



u;2 - 2 — 



1 = fiiA^cs^{f3, m) + ^2B^ds\l3, m) , f3v 
sm{ki)cn{f3,m) = sin(/c2)dn(/3, m) . 



2 sin(A:i) 
cs(/3, m) 



1 - (1 - m)n2B^ 



1 - (1 - m)n2B^' 



2 cos(/c2) 
cn(/3, m) ' 



(65) 



Note that since it is a moving periodic pulse solution, it must not only satisfy the periodicity condition 
()22p but it must also satisfy the periodicity condition (j52p . The periodicity conditions (j22p and (|52p imply 
that Un,Vn are periodic solutions provided there exist integers ni,n2, 713,714 such that the condition ([53]) is 
satisfied. 

In the limit m = 1, these three solution (iii), (iv), and (v) reduce to 



Un = ^exp[— i(a;it — kin + 5i)]sech[/3(n — vt + 82)] , 
Vn = B exp[—i{uj2t — k2n + 53)]sech[/3(n — vt + 82)] , 



(66) 



provided 



ki = k2, sinh2(/3) = (/ii^2 ^ ^2^2) 
wi = 2[1 — cos(/i;i) cosh(/3)] , a;2 = 2 



/i2 

^1 



cos(/ci) cosh(/3) 



, 77/3 = 2sin(/!:i)sinh(/3) . (67) 



Thus the periodic pulse solutions exist provided at least one out of ^2 is positive. For an entirely 
different coupled AL model, solution ([66]) has also been obtained in [33] . 
(vi) Finally, it also admits a coupled periodic kink solution 

Un = ^exp[— 7(tJit — A;i77 + (5i)]-v/msn[/3(77 — vt + 82), m] , 
Vn = B eyip[—i{uj2t — k2n + 83)]^/msn[/3{n — vt + 82), m] , 
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provided 

ki = k2 , wi = 2[1 — cos(A:i)cn(/3, m)dn(/3, m)] , 0J2 = 2 — — cos(A;i)cn(/3, m)dn(/3, m) 
HiA"^ + I22B^ = -sn^{P,m) , /3t; = 2 sin(A:i)sn(/3, m) . (69) 

Thus this solution only exists if at least one out of /Ui, ^2 is negative. Note that since it is a moving periodic 
kink solution, it must not only satisfy the periodicity condition ()25p but it must also satisfy the periodicity 
condition ([52]) . The periodicity conditions (p5]) and ([52|) imply that are periodic solutions provided 

there exist integers ni,n2,ns,n4 such that the condition (|56p is satisfied. 
In the limit m = 1, this solution reduces to 

Un = Aexp[—i{ujit — kin + 5i)] tanh[/3(n — vt + 82)] , 

Vn = B exp[—i{uj2t — k2n + iJs)] tanh[/3(n — vt + 82)] , (70) 

provided 

ki = k2 , oJi = 2[1 — cos(A;i)sech^(/3)] , a;2 = 2 — — cos(A;i)sech^(/3) 
HiA"^ + ^2^^ = - tanh2(/3) , /3v = 2 sm{ki) tanh(/3) . (71) 

While obtaining these solutions, we have made use of several identities for the Jacobi elliptic functions 



5 Solutions of the Coupled Saturated DNLS equations 

We show that unlike the coupled Salerno and the coupled AL case, the coupled saturated DNLS Eqs. ([5]) 
and dni) while they admit Lame polynomial solutions of order one, they do not admit general solutions 
characterized by Eqs. (fT6]) and pT|) or Eqs. (j39]) and (pUj) . In particular, this model does not admit Lame 
polynomial solutions of order two. It is worth noting here that the uncoupled saturated DNLSE model 
does admit Lame polynomial solutions of order one }35j . 

It is easy to check that the coupled Eqs. ([5|) and ([6]) have the following exact solutions in terms of 
Lame polynomials of order one. 
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(i) It admits a coupled mixed pulse-kink solution 



Un = Aexp[-i{ojit + (5i)]dn[/3(n + 62), m] , 

Vn = B exp[-i{uj2t + 63)]^/rnsn[|3{n + 62), m] , (72) 



provided 



oJi = 2 - — , 0J2 = 2 - — , fi2 = /iicn(/3, m) , 

m^^ = , , -1, f^2B' = ^^4%^-^- (73) 



2/iidn(/3, m) ' 2/iidn(/3,m 

For the solution ([73]) . Un,Vn satisfy the boundary condition ([22]) . 
(ii) Another coupled mixed pulse-kink solution is 



tin = ^exp[— i(a;it -|- 5i)]-v/mcn[/3(n -|- 82), m] , 

Vn = B ex.p[-i{uj2t + 63)]^/msn[|3{n + 62), m] , (74) 



provided 



wi = 2 ^, u;2 = 2 , fi2 = /iidn(,9, m) , 

/i2 '^i 1 n2 J^idn2(/3,m) 

2//icn(p,mj 2;Uicn(p,mj 

For the solution ([7i|) . m„,u„ satisfy the boundary condition ([25]) . 

In the limit m = 1, these two solutions (|72p and ()74p go over to the mixed hyperbolic pulse-kink solution 

= Aex.p[—i{u}it + 6i)]sech[(3{n + 62)] , 
Vn = B exp[-i{uj2t + 53)] tanh[/3(n ^2)] , (76) 



provided 



oJi = 2 , UJ2 = 2 , /i2 = ^isech(/3) , 

= = — ^-1. (77) 

2/ii 2/xi cosh(p) 
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(iii) This model also admits two coupled pulse solutions. One solution is 

Un = Aex.p[—i{ujit + 5i)]dn[/3(n + 62), m] , 
Vn = B exp[-i{u}2t + (53)]dn[/3(n + 82), m] , 



provided 



OJl = UJ2 = ^ 



1 



dn(/3, m) 



cn^{/3, m) 
sn^(/3, m) 



cn2(/3, m) 

For the solution (j78|) . Un,Vn satisfy the boundary condition ([3T]) . 
(iv) Another pulse solution is 

Un = A exp[— i(a;it + 6i)]\'^cn[/3{n + 82), m] , 
Vn = B exp[-i(cj2i + (53)]\/mcn[/3(n + (52), m] , 



provided 



U)i = (jJ2 = 2 



fiiiA^ + B^ 



1 



cn(/3, m) 
dn^(/3, m) 
sn^(/3, m) 
dn^(/3, m) 



1^1 = , Ml = 



(78) 



(79) 



(80) 



(81) 



For the solution ()80p . Un,Vn satisfy the boundary condition (|25p. 

In the limit m = 1, these two solutions (j78p . ()80p reduce to the coupled hyperbolic pulse solution 

Un = Aexp[-i{ujit + (5i)]sech[/3(n + ^2)] , 
f„ = B exp[-i{uj2t + 53)]sech[/3(n + 62)] , 



(82) 



provided 

wi = a;2 = -2[cosh(/?) - 1] , z/i = 1^2 , Hi = ^2 
Hi{A^ + B^) = smh\P) . 

Notice that the coupled pulse solutions are admissible only if /ii = /U2 > 0. 



(83) 
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(v) Finally, these coupled equations also admit a mixed kink solution 

Un = Aexp[—i{ujit + 5i)]\/msn[/3{n + 82), m] , 

Vn = B ex.p[-i{LJ2t + 63)]^/rnsn[f3{n + 62), m] , (84) 

provided 

uji = UJ2 = 2[1 — cn(/3, m)dn(/3, m)] , 

ui = U2, Ail = Ai2 , + B^) = -sn2(/3, m) . (85) 

Thus unlike the coupled pulse solution, the coupled kink solution is only valid if /Ui = /i2 < 0. For the 
solution dHU), Un,Vn satisfy the boundary condition (f25|) . 

In the limit m = 1, this solution reduces to the hyperbolic kink solution 

Un = Aeyi.]i[—i{(jJit + Si)] tanh[/3(n + 82)] , 

Vn = B exp[-i(a;2t + h)] tanh[/3(n + 82)] , (86) 

provided 

cji =W2 = 2tanh2(/3), m = fi2 , 1^1 = 1^2, Hi{A^ + B^) = - tanli^{f3) . (87) 

Before ending this section, it might be worthwhile explaining why this model (unlike coupled Salerno 
or coupled AL models) does not admit Lame polynomial solutions of order two. If we look at the field 
equations which follow by using the general PB structure given by Eq. (HD, then it is easily seen that the 
model admits Lame polynomials of order two as solutions provided 

1 + Ai|n„|2 + A2|w„|^ = 0. (88) 

It is easily checked that while this condition can be readily satisfied in both coupled Salerno and coupled 
AL models (where Ai = f^i,X2 = 1^2), this condition can never be satisfied in the coupled DNLS case 
since in that case Ai = A2 = 0. Note, however, that in view of the nontrivial identities for Jacobi elliptic 
functions |34i the coupled DNLS model still admits Lame polynomial solutions of order one. 
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6 Solutions for a Coupled Discrete 0^ Model 

We start from the same continuum coupled 4>^ model for which recently we have obtained Lame polynomial 
solutions of order two [1]. We now show that if we consider the following discrete variant of the same 
model, then it has solutions not only in terms of Lame polynomials of order one but even in terms of Lame 
polynomials of order two, even though the Lame polynomials of order two are not the solutions of the 
uncoupled discrete (j)^ model. 

The field equations for the static coupled continuum model, which we had considered recently [1], are 
given by (modulo a factor of 2 in the definitions of ci, C2, e, /) 

^ = ai(/> - hi(t>^ + 2ci(/)5 + + 2e</.3'02 + 2fH^ , (89) 

^ = aaV- - b2i^^ + 2c2^p^ + d^jcj)'^ + ec/)^ + '^f^P^i'^ ■ (90) 
Let us consider the following coupled discrete model 

1 , , 



n+l + (pn-1 - 20„) = ai(j)n - + dlpl(j)n + [ci(j)'^ + e(/)^Vn + /V'n] [4'n+l + (pn-l] , (91) 



1 e 

■J^{lpn+l + - 2^„) = a2V'n - b2lpn + d(t>liJn + [C2'4'n + ^^'i + f (fii'4''i]Vl'n+l + ^n-l] , (92) 



which in the continuum limit goes over to Eqs. (|89j) and (|90|) . Here h denotes the discreteness parameter. 
Solutions of the Uncoupled Model 

Let us first note that the uncoupled field Eq. (|9ip for field (f) (similar conclusion is also valid for the 
field given by 



1 



Wl + 4>n-l - 2(t)n) = ai0n - bi4>\ + Ci0^[(/)„+i + , (93) 



has three solutions in terms of Lame polynomials of order one. However, Lame polynomials of order two 
do not satisfy the uncoupled Eq. (|93|) . In particular, it is easily shown that 

(t)n = Adn{j3{n + XQ),m], (94) 

is an exact solution to the field Eq. ([93|) provided 



2 

^^/i^cics^(/3, m) = 1, a-i = To 



dn(/3, m) 
cn2(/3, m) 



1 



b{ dn^(/3,m) 



2aiCi cn2(/3, m)[dn(/3, m) — cn2(/3, m)] 
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(95) 



For the solution (j94p . (j)^ satisfies the boundary condition 

Yet another solution to the field Eq. is given by 

4>n = Ay/mcn[l3{n + xo),m\ , 

provided 



2 

cids'^ {(3 , m) = 1, 01 = 777 



cn(/3, m) 
dn^(/3, m) 



5f m^cn^(/3, m) 



2aici dn''(/3,m)[cn(/3,m) -dn''(/3,m)] 



For the solution (I97p . satisfies the boundary condition 



(?J 4i£(m). = <Pn • 

n-i- ^ 

In the limit m = 1, both these solutions go over to the pulse solution 

(pn = Asech[^(n + xq)] , 

provided 

h^A'ci = sinh4(/3) , ai = ^[cosh(/3) - 1] > , - 



/i^ ' 2aiCi cosh(/3) — 1 

The third periodic solution to the field Eq. ()93p is given by 

(j)n = A^/m.sn[f3{n + xo),m] , 

provided 

.4,9 X 2 , ,^ ,^ , , 6? m^cn^(/3, m)dn^(/3, m) 

A^/i^cins^(/3,m) = 1, oi = [cn(/3, m)dn(/3, m) - 1] < , ^ - 



/i^ ' ' ' 2|ai|ci 1 — cn(/3, m)dn(/3, m 

For the solution (I102p . (pn satisfies the boundary condition ([99]) . 
In the limit m = 1, this solution goes over to the kink solution 

4>n = Atanh[P{n + xq)] , 

provided 

/i^^^ci = tanh^(/3) , ai = -^tanh2(/3) < ^1 _ 4 



/i2 ' 2|ai|ci sinh2(2/3) 
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Let us now discuss the solutions of the coupled Eqs. (j9T|) and (|92]) . It turns out that as in the coupled 
AL case, the (f)^ coupled equations have solutions satisfying the ansatz similar to (|17|) (but no solutions 
satisfying the ansatz similar to (|4Up ). and also solutions in terms of Lame polynomials of order one (by 
making use of the identities for the Jacobi elliptic functions [34]). 
Solutions of the Coupled Model Satisfying Ansatz Similar to ()17|) 

On substituting the ansatz 

4>l + ai)l = h, a,b > 0, (106) 

(which is similar to the ansatz (fT7|) l in the coupled field Eqs. (pTj) and (j92]) . we find that such solutions 
exist provided 

ci = C2 = / = - , bi = b2 = -d, ai = a2, ci/i^6^ = 1 , 

a = l, -1 + ^ = ^^- (107) 

This is a rather general ansatz and there are several solutions of this type which exist for this model. 
Lame polynomial solutions of order one 

(i) One solution is 

(j)n = Adn[l3{n + C2),m] , ^lJn = B^/msn[(3{n + C2),rn] , (108) 
provided Eq. ()107p is satisfied and further 

b = A^, = 3"^ , ai- biA^ + 2ci = . (109) 

Note that the width (3 is completely arbitrary. For this solution, (pn,ipn satisfy the boundary condition 
given by Eq. ([2^ with (j)n,ipn replacing fn,gn respectively. 

(ii) Another solution is 

4>n = A^/mcn[f3{n + C2),'m] , ipn = B^/msn[|3{n + C2),'m] , (110) 
provided Eq. (|107p is satisfied and further 

b = mA^ , A^ = B^ , ai - mbxJ^ + Im^cxA^ = . (Ill) 
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For this solution, (l)n,ipn satisfy the boundary condition given by Eq. (j25|) with (l)n,'>Pn replacing fn,gn 
respectively. 

In the limit m = 1, both these solutions go over to the hyperbolic solution 

/„ = Asedi[f3{n + cs)] , 9n = B tanh[/3(n + cs)] , (112) 

Lame polynomial solutions of order two 

(iii) One solution is given by 

(pn = Adn'^[/3{n + C2),m] + B , Vn = -^V^sn[/3(n + C2), m]dn[/3(n + C2), m] , (113) 

provided Eq. (1107P is satisfied and further 

42 

6= — , A^ = F'^, A = -2B, ai - 61^2 + 2ci 5^ = 0. (114) 

For this solution, (pnii^n satisfy the boundary condition given by Eq. (j22p with (j)n,il^n replacing fn,gn 
respectively. 

(iv) Another solution is 

4>n = Adn'^[(3{n + C2),m] + B , = -?^"isn[/3(n + C2), m]cn[/3(n + C2), m] , (115) 
provided Eq. (I107p is satisfied and further 

^^m_^ ^2^i?2^ {2-m)A = -2B , 8ai - 261/71^^2 + cim^F^ = . (116) 

For this solution, <j)n,tpn satisfy the boundary condition given by Eq. (fSTI) with <j)n,tpn replacing fn,gn 
respectively. 

In the limit m = 1, both the solutions (jll3p and (jllSp . go over to the hyperbolic solution 

= ^sech2[/3(n + C2)] + 5, = i^tanh[/3(n + C2)]sech[/3(n + C2)] . (117) 

(v) Apart from these, several other solutions are possible. For example one can have the following 
nonperiodic solution 

^ , Bn , , 

i;n = r. o ' (118) 
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provided Eq. ()107p is satisfied and further 

6 = ^2, A'^ = B^. (119) 

(vi) Yet another solution is 

^n = Acos[l3{n + C2)], = ^sin[/3(n + C2)] , (120) 
provided Eq. ()107p is satisfied and further 

b = A'^, A^ = B^. (121) 

For this solution, (pn,'4'n satisfy the boundary condition given by Eq. (|55|) with (pn,'4'n replacing fn,gn 
respectively. 

Solutions following from identities for Jacobi elliptic functions 

Using the identities for the Jacobi elliptic functions [34J, we now show that there are six Lame polynomial 
solutions of order one to the coupled field Eqs. ([9T]) and ([92]) . 
Solution 1: It is not difficult to show that 

(pn = Adn[l3{n + Xo),m] , ipn = B^/msn[|3{n + xo),m] , (122) 

is an exact solution to the coupled field Eqs. (i9T]l and (f92]) provided 

61^2 + dB^ = 2(ci + fB^ - eA'^B^)ds{P, m)ns(/3, m) , (123) 

ai + ^ + dB^ + 2{eA^B'^ - 2/5^)ds(/3, m)ns(^, m) = 2(ciA^ + /fi^ - eA^B'^)ds(l3, m)ns(/3, m)cs^{l3, m) , 

(124) 

^-fB^ + {eA^B'^ - 2fB^)cs^{p, m) = {ciA^ + - eA'^B^)cs\p, m) , (125) 
62^^ + dA^ = -2 (^C2B^ + ^A^ - 2fA^B^^ ds(/3, m)cs(/3, m) , (126) 

a2 + ^ + dA'^- 2(2fA^B^ - eA^)ds(/3, m)cs(/3, m) = 2(02^^ + -A^ - 2fA'^B^)ds((3, m)cs(/3, m)ns2(/3, m) , 
ft^ 2 

(127) 

_ £^4 _ (2/^2^2 - eA^)ns\(3, m) = ( C2B^ + -^^ - 2fA^BA ns^(/3, m) . (128) 
/i^ 2 V 2 / 
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For this solution, (l)n,ipn satisfy the boundary condition given by Eq. (j22j) with (l)n,'>Pn replacing fn,gn 
respectively. 

These equations are of course trivially satisfied if relations (|107p and (|lU9p are satisfied. It is worth 
pointing out that while relations (|lU7p and (|lU9p are sufficient so that (|122p constitutes an exact solution 
to the coupled Eqs. (|9ip and (|92p . it is not obvious if relations (|lU7p and (|lU9p are also necessary. The 
necessary relations are as given by Eqs. (|12c{p to (|128p . 
Solution 2: Another solution is given by 

(j)n = A^/rncn[|3{n + Xo),rn] , ipn = B^/msn[f3{n + xo),m] , (129) 

which is an exact solution to the coupled field Eqs. ([9T]) and (j92]) provided 

biA^ + dB'^ = 2{cxA^ + - eJ^B'^)z&{^P, m)ns(/3, m) , (130) 

ax^^^mdB'^^2m{eA^B'^-1fB^)c^{^,m)u^{^,m) = 2(01^^ + fB^-e A^B'^)csi (3, m)iisip,m)ds^ {13, m) , 

(131) 

^ - m^fB^ + mieA^B^ - 2fB^)ds\p, m) = (ci^^ + /^^ - eA^B^)ds\p, m) , (132) 
b2B^ + dA^ = -2 (c2B^ + ^A"^ - 2fA^BA ds(/3, m)cs(/3, m) , (133) 



a2+^+mdA^-2i2fA^B^-eA^)cs{P,m)csil3,m) = 2 (^ca^^ + - 2fA^B'^^ ds(/3, m)cs(/3, m)ns2(/3, m) , 

(134) 



^ - -m'A'' - {2fA'B' - eA^)ns\t3, m) = ca^^ + -A"" - 2fA'B' ns\f3, m) . (135) 
2 \ 2 J 

For this solution, (j)n-,'4'n satisfy the boundary condition given by Eq. ([251) with (pn,'4'n replacing fn,9n 

respectively. 

These equations are trivially satisfied if relations (jlUZp and (jllip are satisfied. 

In the limit m = 1, both the solutions (I122p and ()129p go over to the hyperbolic soliton solution 

4>n = Asech[(3{n + xq)] , ^„ = i? tanh[/3(n + xq)] , (136) 

provided relations (|123p to (jl28p with m = 1 are satisfied. 
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Solution 3: It is not difficult to show that 

(f>n = Adn[(3{n + Xo),m] , ^pn = B^/mcn[|3{n + xo),rn] , (137) 
is an exact solution to the coupled field Eqs. ([9T]) and (j92]) provided 

biA^ - dB^ = 2(ci^^ + /S^ + eA^B^)ds{P, m)ns(/3, m) , (138) 
2 

«i + T9 - (1 - m)dB'^ - 2(1 - m)(eyl^B2 + 2/S^)ds(/3, 

= 2(ci^^ + /B^ + eA^B'^)ds{P, m)ns(/3, m)cs2(/3, m) , (139) 



- (1 - m)2/54 - (1 - m)(e^2^2 ^ 2fB'^)cs\p, m) = (ci^^ + /^^ + e^252)pg4(^^ ^ ^^40) 

62^^ - dA^ = -2 (^C2B^ + |yl^ + 2fA^B^^ ns(/3, m)cs(/3, m) , (141) 

02 + ^ + (1 - m)(i^2 + 2(1 - m){2fA^B^ + eA^)ns(/3, m)cs(/3, m) 

= 2 (^C2B^ + + '^f^^B^^ ns(/3, m)cs(/3, m)ds2 (;5, m) , (142) 

^ - |(1 - m)2^^ + (1 - m)i2fA^B^ + eA^)ds\l3, m) = (ca^^ + + 2fA^B^)As\p, m) . (143) 



For this solution, 4>mipn satisfy the boundary condition given by Eq. (j22|) with (pmipn replacing fn,gv 
respectively. 

Solution 4: Another solution to the coupled Eqs. ([9T]) and ([92|) is given by 



provided 



(j)n = Adn[/3{n + Xo),m] , ipn = Bdn[f3{n + xo),m] , (144) 

biA^ - dB^ = 2(ci^^ + /B^ + eA^B^)ds{P, m)ns(/3, m) , (145) 

ai + ^ = 2iciA^ + fB^ + e^^^^^^jsp^ m)ns(/3, m)cs2(/3, m) , (146) 

-3^ = {ciA^ + fB^ + eA^B^)cs\l3, m) , (147) 

62^^ - = 2 (^C2S^ + + 2fA^B'^^ ds(/3, m)ns(/3, m) , (148) 

a2 + = 2 (^C2B^ + + 2fA^B^^ ns(/3, m)ds(/3, m)cs'^{f3, m) , (149) 
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^ = (^C2B^ + ^ 2/^2^2^ cs\p, m) . (150) 



For this solution, 4'n,'>Pn satisfy the boundary condition given by Eq. (j3T|) with (pmipn replacing fn,9v 
respectively. 

Solution 5: Yet another solution to the coupled Eqs. ([9T]) and (j92]) is given by 



provided 



(pn = Ay/rncn[l3{n + Xo),m] , ipn = B\/mcn[/3{n + xo),m] , (151) 



biA^ - dB^ = 2(ci^^ + + eA^B^)cs{(3, m)ns(/3, m) , (152) 
2 

ai + — = 2(ci^^ + /B^ + eA^B^)cs{/3, m)ns(/3, m)ds2(/3, m) , (153) 

^ = (ci^^ + fB^ + eA^B^)ds\(3, m) , (154) 

b2B^ -dA^ = 2 (^C2B^ + ^A^ + 2fA^B^^ cs(/3, m)ns(/3, m) , (155) 

02 + = 2 (^C25^ + + 2/^2^2^ ns(/3, m)cs(/3, m)ds2(/3, m) , (156) 



^2 - yC25* + -A"" + 2/^^5^ j ds*(/3, m) . (157) 

For this solution, cpm'^'n satisfy the boundary condition given by Eq. ([25]) with cpuji^n replacing fn,gn 
respectively. 

In the limit m = 1, all three solutions given by (I137p . (11440 and (I15ip go over to the hyperbolic soliton 
solution 

4>n = j4sech[/3(n + xq)] , ipn = Bsech.[(3{n + xq)] , (158) 

provided 

2 2 

ai = 02 = ^[cosh(/3) - 1] > , ftivl^ - dB'^ = 62^^ - = sinh2(/3) cosh(/3) , (159) 

ciA^ + eA^B^ + fB' = c^B^ + 2 /A^i?^ + ^ ^mh^ ^ 

2 /i^ 

Solution 6: Finally, another solution to the coupled Eqs. (|9ip and (j92p is given by 

= A-v/msn[/3(n + Xo), m] , = i3-v/msn[/3(n + xq), m] , (161) 
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provided 



biA^ - dB"^ = 2(ci^^ + fB^ + eA^B'^)cs{P, m)ds(/3, m) , (162) 

ai + = 2{ciA^ + fB^ + e^2^2vg(^^ m)ds(/3, m)ns2(/3, m) , (163) 

62^2 - ^ 2(c25^ + + 2/^2^2)cs(/3, m)ds(/3, m) , (165) 

02 + = 2(c25^ + + 2fA^B'^)ds{l3, m)cs(/3, m)ns2(/3, m) , (166) 

^ = (C2i?^ + + 2/A2s2)ns4(/3, ,„) . (167) 



For this solution, <?!)n,V'n satisfy the boundary condition given by Eq. ([251) with 4'n,'>Pn replacing fn,gn 
respectively. 

In the limit m = 1, this solution goes over to the hyperbolic soliton solution 

4>n = A tanh[/3(n + xq)] , i^n = B tanh[/3(n + xq)] , (168) 

provided 

ai = 02 = -To [tanh^ (/?)]< 0, 61^" - = b2B^ - dA^ = . ^ , (169) 

n,^ cosh (p) 

ci + eyl2i?2 + ^ ^^^4 ^ 2 fA^B^ + -yl^ = ^^"^'^^^ . (170) 

2 n,^ 

7 Solutions for a Coupled Discrete cj)^ Model 

We start from the same coupled static discrete field equations as in our recent paper [29] for which we 
had obtained six solutions in terms of Lame polynomials of order one. We now show that the same model 
also admits two Lame polynomial solutions of order two, even though they are not the solutions of the 
corresponding uncoupled problem. Let us start from the field equations considered in [29J 

^{(t>n+l + 4>n-l - 2</)„) - 2ai(t>n - [2^1 + 7V'n]['/'n+l + (l>n-l] = , (171) 
^(V-n+l + i'n-l - 2^n) " 2a2V'n " [2/32V'n + 7</'^][V'n+l + i'n-l] = . (172) 
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Let us now discuss the solutions of the coupled Eqs. (jl7ip and (|172p . It turns out that as in the 
coupled (j)^ case, the (j)'^ coupled equations have solutions satisfying ansatz similar to the one given by Eq. 
p7|) (but no solutions satisfying ansatz similar to (|4U|) ). Further, they also have solutions in terms of Lame 
polynomials of order one which we have already discussed in Note that these solutions were obtained 
by making use of the identities for the Jacobi elliptic functions [34| . 
Solutions satisfying ansatz similar to ()17l) 

On substituting the ansatz as given by Eq. (|1U6|) (which is similar to the ansatz given by Eq. (|17|) ) in 
the coupled field Eqs. (|171|) and (|172|) . we find that such solutions exist provided 

2/3i = 2/32 = 7, ai = a2 = -^, a = l. (173) 

This is a rather general ansatz and there are several solutions of this type which exist in this model. 
Lame polynomial solutions of order one 

(i) One solution is 

(/>n = ^dn[/3(n + C2),m] , t/;^ = -B\/msn[/3(n + C2), m] , (174) 
provided Eq. ()173p is satisfied and further 

Note that the width /? is completely arbitrary. For this solution, (pn,i^n satisfy the boundary condition 
given by Eq. (p2]) with (j)n,ipn replacing fn,gn respectively. 

(ii) Another solution is 

(j)n = A^/mcn[|3{n + C2),rn] , tpn = B^/msn[|3{n + C2),'m] , (176) 
provided Eq. (|173p is satisfied and further 



For this solution, 4'n,'>Pn satisfy the boundary condition given by Eq. (j25|) with 4'n,'>Pn replacing fn,9v 
respectively. 
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In the limit m = 1, both these solutions go over to the hyperbolic solution 

(j)n = Asech[/3{n + C2)] , ^pn = B tanh[/3(n + C2)] , (178) 

Lame polynomial solutions of order two 

(iii) One solution is given by 

4)n = Adn'^[(3{n + C2),m]+B, = i^\/^sn[/3(n + C2), m]dn[/3(n + C2), m] , (179) 
provided Eq. ()173p is satisfied and further 



For this solution, (t)n-,i^n satisfy the boundary condition given by Eq. (j25|) with 4'n,'>Pn replacing fn,9n 
respectively. 

(iv) Another solution is 

= j4dn^[/3(n + C2), m] + B , = -?^"T-sn[/3(n + C2), m]cn[/3(n + C2), m] , (181) 
provided Eq. (1173P is satisfied and further 

2*2 

^^m_^ A2^i?2^ {2-m)A = -2B. (182) 

For this solution, (pnjfpn satisfy the boundary condition given by Eq. (|3T]) with (pnji^n replacing fn,gn 
respectively. 

In the limit m = 1, both solutions (|179p and (|18ip go over to the hyperbolic solution 

(/)n = ^sech2[/3(n + C2)] + -B, ^„ = Ftanh[/3(n + C2)]sech[/3(n + C2)] . (183) 

(v) Apart from these, several other solutions are possible. For example one can have a nonperiodic 
solution 

A , Bn , , 

Vn = r. o ' (184) 



provided Eq. (|173p is satisfied and further 



b = A'^, A'^ = B^. (185) 
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(vi) Yet another solution is 



= Acos[l3{n + C2)] , tpn = Bs'm[/3{n + C2)] 



(186) 



provided Eq. (I173P is satisfied and further 



b = A\ A^ = B\ 



(187) 



For this solution, (pn,'4'n satisfy the boundary condition given by Eq. ([3T]) with (j)n,'^n replacing fn,gn 
respectively. 

8 Summary 

In this paper we have shown that for a number of coupled discrete models, e.g. coupled Salerno, coupled 
Ablowitz-Ladik, coupled saturated nonlinear Schrodinger equation, coupled (j)^, coupled (p^, while the 
uncoupled equations do not admit solutions in terms of Lame polynomials of order two, the coupled 
models do admit such solutions. These solutions (with appropriate boundary conditions) have relevance 
in physical contexts ranging from ferroelectric [7, 8, 9] to multiferroic [4, 5, 6] materials to the models in 
field theory [10, 14] as well as for various discrete contexts [11, 12, 13]. 

The stability of various solutions found here remains an open issue to be explored numerically, par- 
ticularly some solutions have an arbitrary soliton width. In addition, the scattering of solitons of various 
discrete models is an important issue with these static solutions boosted with a certain velocity. Similarly, 
the Peierls-Nabarro (discreteness) barrier for the solutions remains to be explored. Given the solutions in 
terms of Lame functions of order one and two, it is then worth enquiring if one considers coupling of three 
discrete fields, would they admit solutions in terms of Lame polynomials of order three? And if true, can 
one generalize it to the case of coupled fields? We hope to address these issues in the near future. 
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